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摘要 : JERI ESFE (nonlinear Schrödinger equation, NLSE) 是 多 类 物理 现象 所 遵守 规律 
的 共同 数学 形式 。 在 光学 系统 中 ，NLSE 可 以 描述 激光 脉冲 在 介质 中 传播 的 时 空 演化 ， 在 天 文学 
动 暗物质 纤维 等 天 体 和 观测 现象 (BS EXE TS — HIA TT FE AR Zt, SPE), 
Klein-Gordon 系 统 , EKG)。 从 仿真 原理 、 


中 ， 可 以 描述 引力 透镜 、 水 


EXP, modu, 


甚至 可 以 巧妙 地 用 来 类 比 强 引力 场 现象 (所 谓 爱 
光路 设计 、 实 验 结果 、 天 文 诠释 等 方面 ， 详 细 介绍 了 非 线性 光学 仿真 实验 在 天 文学 领域 的 应 用 和 


ERY, KRE, 
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FREN, HeH 
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进展 ， 并 简要 展望 了 未 来 
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1 前 言 


类 比 引力 ”是 一 种 类 比 研究 方法 ， 通 过 料 
态 甚 至 量子 信息 等 ) 实验 来 研 
现象 。 类 比方 法 的 坚实 基础 在 于 : 引力 现象 和 仿真 实验 由 
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IE 


DES 


(SEAR HERE SITE (NLSE)”， 既 可 以 在 光学 中 描述 激光 脉冲 在 介质 中 传 所 


演化 ， 也 可 以 描述 引力 透镜 、 波 动 暗物质 (wave dark matter, waveDM)” 等 天 体 和 天 文 观 
测 现象 ， 甚 至 可 以 巧妙 地 用 来 类 比 强 引力 场 E 
要 有 如 下 两 方面 : (1) 通过 高 精度 的 、 可 调控 的 实验 ， 人 们 可 精 引 


以 霍金 辐射 这 种 黑洞 视界 阶 
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观测 数据 。(2) 希望 可 以 为 天 体 物 天 


f 究 或 者 说 仿真 难以 观测 、 久 


展 非 线 性 光学 实验 ， 仿 真 各 种 波动 暗物质 结构 的 演化 之 前 景 。 
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心 设计 的 非 引力 系统 (如 流体 、 光 学 、 凝 聚 
以 实验 的 天 体 物理 学 和 宇宙 学 引力 


同样 的 数学 方程 来 描述 。 例 如 本 文 
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仿真 引力 的 实验 ， 意 义 主 
旧地 探测 和 研究 强 引 力 场 下 
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这 样 的 效应 ， 特 别 是 难以 获得 一 些 高 质量 观测 数据 的 天 文 现象 。 
近 的 量子 效应 为 例 ， 直 接 探测 难度 很 大 ， 
LE 学 和 宇宙 学 中 的 研究 对 象 提供 启发 性 的 新 思路 。 这 是 


至 今 还 没有 令 人 信服 的 
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为 天 文学 系统 内 在 自由 度 往往 很 多 ， 而 且 描述 这 些 系统 的 数学 方程 高 度 非 线性 ， 所 以 不 可 能 
解析 求解 ， 甚 至 很 难 数值 求解 。 因 此 我 们 希望 在 实验 室 做 仿真 实验 来 理解 天 文学 系统 ;特别 
是 ， 由 于 天 文学 对 象 的 演化 时 标 很 长 ， 人 类 不 可 能 跟踪 观测 它们 的 演化 过 程 ， 因 此 来 自 仿真 
实验 的 认识 就 尤为 重要 ( 详 见 4.2 节 和 第 5 章 )。 当 然 ， 对 于 同一 个 艰深 复杂 的 数学 形式 ， 从 
不 同 的 角度 来 理解 和 处 理 ， 往 往 在 具体 的 不 同方 面 带 来 独特 的 便利 性 ， 因 此 引力 理论 也 可 以 
帮助 这 种 在 实验 室 开展 的 精细 实验 ， 即 类 比 是 双向 的 (如 用 引力 领域 的 虫 洞 来 理解 复杂 的 量 
子 纠缠 ， 指 导 实 验 设计 ”)。 

自从 广义 相对 论 发 明 以 来 ， 人 们 就 提出 了 一 些 理论 思路 和 模型 来 理解 弯曲 时 空 ， 如 基于 
光线 在 可 变 折线 率 介质 传播 的 几何 光学 方程 (所 谓 Gordon 度 规 ”)。 一 般 认 为 ， 现 代 类 比 引 
力 的 标志 性 突破 是 Unruh” 在 1981 年 提出 的 “声学 黑洞 模型 >， 用 类 比 的 方式 来 “观察 ” 霍 
4. Unruh 用 流体 流动 中 的 声波 运动 ， 来 类 比 强 引 力 场 中 无 质量 粒子 的 运动 (即刻 画 出 
弯曲 时 空 的 测 地 线 )。 他 具体 分 析 了 这 两 种 现象 所 具有 的 相同 的 数学 形式 ， 如 在 跨 声 速水 流 
中 的 声波 视界 对 应 的 就 是 黑洞 的 视界 ， 声 波 视界 附近 发 出 的 热 谱 形式 的 声波 对 应 于 黑洞 的 
霍金 辐射 。 在 随后 的 几 十 年 中 ， 提 出 的 有 关 仿 真实 验 的 理论 方案 层出不穷 ， 设 想 的 实验 手段 
从 普通 流体 延伸 到 光学 系统 、 玻 色 - 爱 因 斯 坦 凝 聚 体 (Bose-Einstein condensate, BEC) 这 样 
的 量子 系统 。 但 真正 开展 的 仿真 实验 ， 一 直 等 到 21 世纪 10 年 代 末期 才 出 现 。 

类 比 引力 的 仿真 实验 最 先是 在 光纤 (以 及 水 槽 ) 中 开展 ””， 随 后 两 年 又 有 了 BEC X 
验 ”， 一 开始 是 瞄准 最 令 人 感 兴 趣 的 目标 (基本 上 在 2017 年 前 一 直 如 此 ): 黑洞 的 霍金 辐 
射 ”™”。 后 来 很 自然 地 就 扩大 了 类 比 引力 实验 的 对 象 ， 如 仿真 旋转 黑洞 的 性 质 ， 以 及 提取 黑 
洞 旋 转 能 量 的 彭 罗斯 过 程 等 各 种 天 体 物理 学 和 宇宙 学 现象 “~”。 详 细 的 历史 发 展 背 景 介绍 ， 
请 参见 文献 四 B, 22, 到 。 最 近 几 年 对 于 弱 引 力 场 的 天 体 和 现象 (例如 引力 透镜 、 潮 汐 效 应 ， 
以 及 非 相对 论 性 的 类 似 玻 色 星 天 体 ) 也 开展 了 光学 仿真 实验 ， 特 别 是 对 于 天 文系 统 长 时 标 演 
化 过 程 开展 了 仿真 实验 。 在 天 文 传统 中 ， 这 类 对 象 正 是 普 裔 而 重要 的 天 文 观测 对 象 。 正 如 第 
一 段 所 述 ， 天 体 演化 过 程 的 仿真 实验 对 天 文学 家 启示 很 大 。 

本 文 主要 综述 激光 实验 仿真 引力 系统 的 最 新 研究 进展 。 并 且 ， 对 于 当前 天 文学 界 展开 热 
烈 讨论 的 波动 暗物质 ， 我 们 简单 介绍 光学 仿真 实验 的 前 景 。 与 流体 实验 相 比 ， 激 光 实验 具有 
精细 调控 ， 精 度 高 的 优点 ;如 水 槽 实验 中 ， 类 比 霍金 辐射 的 声波 信号 与 水 中 的 各 种 波动 (如 
热 噪声 ) 相 比 太 微弱 。 跟 冷 原 子 实验 相 比 (如 BEC)， 激 光 实 验 的 实现 难度 低 ， 成 本 低 。 

本 文 结构 如 下 : 第 2 章 ， 介 绍 光 学 非 线性 薛 定 廖 方程 (NLSE) 系统 仿真 引力 系统 的 原理 
(仿真 的 技术 思路 ); 第 3 章 和 第 4 章 ， 从 实验 的 光路 设计 、 实 验 结 果 和 所 对 应 的 天 文 诠释 三 
个 方面 ， 分 别 对 仿真 强 引力 系统 、 弱 引力 系统 的 光学 实验 进行 详细 介绍 ; 第 5 章 总 结 全 文 ， 
并 展望 利用 光学 实验 研究 波动 暗物质 (waveDM) 多 种 结构 的 动力 学 演化 过 程 。 


2 光学 NLSE 系统 仿真 引力 系统 的 原理 


9 


z2 


在 光学 系统 中 ， 假 设 光 场 沿 轴 z fed. E = ple, y, zje, 95 y 满足 


« 


42 天 文学 进展 42 卷 


P 即 光 场 包 络 va, y, 2) 沿 轴 z 是 慢 变 时 ， 其 动力 学 方程 可 以 由 NLSE 描述 中， 形式 
a RB. 
OW lea, kAÀn, | 
"8s ak V Vs no pum (1) 
2 2 
其 中 ， = ta 为 二 维 拉 普 拉 斯 算 子 ，m 为 介质 折射 率 的 线性 部 分 ，Anm 取决 于 介 
质 的 非 线性 极 化 强度 ， 波 数 大 二 2xno/ 和 。 从 数学 角度 来 看 ，NLSE 通常 没有 解析 解 ， 而 且 


数值 求解 一 般 也 十 分 复杂 。 近 年 来 非 线 性 光学 飞速 发 展 ， 光 学 实验 成 为 模拟 求解 NLSE 的 
一 种 有 效 方法 。 如 第 1 章 所 介绍 ， 通 过 非 线性 光学 实验 ， 现 已 成 功 地 开展 了 天 文 对 象 的 实 
验 室 仿真 ， 不 但 针对 引力 透镜 、 引 力 红 移 和 潮汐 效应 以 及 旋转 玻 色 星 等 一 系列 重要 的 天 文 观 
测 对 象 ， 而 且 针 对 黑洞 的 霍金 辐射 、 绢 罗斯 超 辐 射 等 难以 观测 的 EKG 现象 。 

为 了 仿真 引力 系统 ， 需 要 将 光学 系统 NLSE 中 的 各 个 参数 与 引力 系统 类 似 方程 的 参数 
一 一 对 应 。 根 据 引力 系统 的 强 弱 ( 即 是 否 为 相对 论 性 )， 仿 真 所 基于 的 具体 对 应 原理 (技术 方 
R) 并 不 相同 ， 下 面 我 们 分 别 进行 介绍 。 
2.1 ” 强 引 力 系 统 的 仿真 原理 

在 激光 实验 仿真 黑洞 这 样 的 强 引 力 系统 时 ， 我 们 利用 Unruh 提出 的 “声学 黑洞 模型 ”， 
将 激光 看 作 是 一 种 光子 流体 ， 用 光子 流体 中 微小 扰动 的 传播 (密度 扰动 和 相位 扰动 ， 以 后 为 
了 叙述 方便 ， 简 称 声波 ) 类 比 弯 曲 时 空中 无 质量 粒子 的 传播 ”。 无 质量 粒子 即 强 引 力 场 中 的 
测试 粒子 ， 它 的 运动 描述 了 弯曲 时 空 的 度 规 。 

在 弯曲 时 空中 ， 描 述 一 般 性 的 标量 粒子 运动 的 方程 为 带 时 空 度 规 的 相对 论 性 (高 速 的 ) 
量子 力学 方程 ， 即 Einstein Klein-Gordon (EKG) 方 程 系统 : 


— 


— [16 -- m?$ — 0 . (2) 


达 朗 贝尔 算 符 口 = "VV i> Hip, g^" 为 时 空 度 规 ，m 为 粒子 质量 。 在 研究 历史 上 ， 这 种 
在 自 引力 ( 爱 因 斯 坦 引力 场 ) 束缚 下 的 、 具 有 相对 性 运动 速度 的 标量 场 ， 称 为 Klein-Gordon 
Geon 四 。 本 文 基本 上 不 涉及 这 种 EKG 系统 。 本 节 中 我 们 主要 借用 m = 0 的 EKG 系统 来 
说 明 引 力度 规 ， 而 对 相对 论 性 标量 场 的 性 质 和 运动 不 多 介绍 ; 在 2.2 节 我 们 只 考虑 EKG 的 
非 相 对 论 极限 (低速 十 牛顿 引力 版 本 )， 即 NLSE4-Poisson 方程 。 

下 面 我 们 简单 介绍 ， 光 场 中 的 声波 扰动 如 何 与 强 引 力 场 中 的 无 质量 标量 粒子 的 运 
动 相对 应 。 把 光学 NLSE 方程 ( 见 式 M ) 变换 到 流体 力学 表象 ， 即 Madelung 变换 : 
y — pe, KB. Jp 和 分 别 为 光 场 包 络 的 振幅 和 相位 。 并 且 假 设 折射 率 与 光 强 的 关系 
为 : An=—71 = wv, Xs y 为 介质 的 非 线 性 系数 。 则 式 四 ) 转化 为 一 系列 流体 力学 


?此 处 公式 中 的 记号 与 文献 回 中 的 相应 公式 (该 文 的 公式 3) 略 有 不 同 ， 对 应 关系 如 下 : ASCH pb 与 该 文 的 A 对 应 ， 本 文 
的 大 等 于 该 文 的 moko， 其 中 mo 为 介质 的 折射 率 ，ko = = 为 光 在 真空 中 的 波 数 。 
0 
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中 的 方程 : 
Bip V (pv) — 0 , (3) 
1, c^ c2 V2p!/2 
| 一 0 . 4 
am 2" ns ^ 2k?nà p? 4 9 


st (B) 和 (E) 分 别 为 流体 力学 的 连续 性 方程 和 欧 拉 方 程 ， 为 光子 体系 的 流体 描述 ， 我 们 以 后 
称 之 为 光子 流体 。 其 中 ，p 对 应 着 流体 的 密度 ， 流体 的 速度 "一 7 Vo = Vy 对 应 流体 
的 速度 势 。 同 经 典 流体 相 比 ， 式 @ 中 最 后 项 为 量子 粒子 (严格 多 说 是 波动 性 ) 独 有 的 正 
力 项 ， 一 般 称 之 为 量子 压 ， 它 来 自 于 光学 方程 中 的 衍射 项 ， 是 不 确定 性 原理 所 导致 的 结果 ， 
当 光子 流体 的 密度 变化 缓慢 时 ， 式 @) 最 后 一 项 量子 压 可 以 忽略 。 

Unruh "在 E 导 经 典 流体 中 声波 运动 的 声学 度 规 时 ， 假 设 流体 是 无 笑 、 正 压 ， 且 流动 是 
无 旋 的 ， 这 三 个 条 件 在 一 般 流体 中 很 难 满足 。 但 在 光子 流体 中 ， 式 (0) 欧 拉 方 程 是 理想 流体 
(无 粘 滞 ) 的 动量 方程 ， 即 光子 流体 本 身 大 致 上 可 当 作 是 “ 超 流体 ”。 而 且 在 我 们 这 种 场景 中 
流体 密度 变化 不 大 ， 量 子 压 可 忽略 ， 那 么 光子 流体 的 压强 只 是 密度 的 函数 ， p= oe 为 
正 压 流体 ， 并 且 在 理想 流体 中 ， 流 动 是 势 流 ， 无 旋 这 一 条 件 也 角 — 
心 处 的 奇 点 )。 因 此 ， 光 子 流体 是 比 普通 流体 更 加 理想 的 实验 平台 。 

考虑 用 线性 微 扰 的 方法 处 理光 子 流体 中 的 微 扰 ， 设 p= poten. +O (2), m = mo tem + 
O(e)。po Fil mg 分 别 为 背景 流 场 的 密度 和 速度 势 。 代 入 式 (B) 和 图 ， 得 到 扰动 (p1, m) 的 
动力 学 方程 : 


» 


Opi +V : (PoV  p1vo) =0 , (5) 


cy 
On + Vm -vo + we 0 . (6) 
0 


将 式 (B) 代入 式 (B) 中 ， 可 得 到 只 关于 六 的 方程 : 


v 
0 机 (Om + vo va) +V. [ovn — : (On + vo: Vn)| 20 , (7) 


S 


其 g = SEI) om, 为 声波 在 流体 中 的 传播 速度 ，wo — Vm， 为 背景 流 场 的 流 
速 。 现 在 引入 等 MUN gu 将 式 四 ) 改写 为 : 


9, (V. —-gg""0,m) —0 ， (8) 


1 
"7 Vg 
度 规 的 协 变形 式 gu, 为 : 
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其 中 ，w 和 ve 分 别 为 流速 的 径 向 和 角 向 分 量 ，vo2 = v2 ++ vo 

由 式 加 和 (B) 可 以 看 出 ， 在 引入 等 效 度 规 gu 后， 声波 在 光子 流体 中 的 运动 方程 形式 
上 与 无 质量 标量 粒子 在 引力 场 中 的 运动 方程 相同 ( 即 式 (B) P m = 0)。 不 过 需要 强调 的 是 ， 
我 们 在 上 述 推导 中 有 两 个 假设 : 首先 ， 式 团 中 最 后 一 项 量子 压 要 能 够 忽略 ， 这 就 要 求 光子 
流体 的 密度 变化 必须 十 分 缓慢 ， 具 体 而 言 ， 光 子 流体 中 激发 的 声波 波长 足够 长 ， 则 声波 造成 
的 密度 变化 可 以 忽略 ; 其 次 ， 介 质 折射 率 的 变化 An = -yl ， 这 就 要 求 非 线性 效应 是 局 半 
的 (光学 领域 称 作 Kerr 效应 )， 并 且 介质 是 自 散 焦 的 ， 即 介质 的 折射 率 与 光 强 成 反比 。 只 
上 述 两 个 假设 都 满足 时 ， 才 能 在 光学 实验 中 模拟 黑洞 这 样 的 强 引 力 系统 。 

2.2 985|7] SPE 系统 的 仿真 原理 

在 2.1 节 中 ， 我 们 对 仿真 强 引 力 场 的 思路 (技术 方案 ) 进行 了 详细 的 介绍 。 在 本 节 ， 我 
们 介绍 弱 引 力 场 系统 (牛顿 引力 ， 或 者 说 平 直 时 空 )， 而 且 只 考虑 牛顿 自 引 力 下 的 、 粒 子 运动 
速度 远 低 于 光速 的 标量 场 系统 。 这 种 牛顿 引力 的 低速 标量 场 ， 近 年 来 在 天 文学 中 引起 了 较 多 
的 关注 ， 例 如 波动 暗物质 理论 ”， 我 们 在 第 4 章 和 第 5 章 中 具体 介绍 。 

如 在 2.1 节 所 介绍 ， 这 种 系统 是 EKG 的 非 相 对 论 极限 (低速 + 牛顿 引力 版 本 )。 在 牛顿 
引力 情形 下 ， 爱 因 斯 坦 场 方程 退化 为 泊 松 方程 。 在 粒子 运动 速度 远 低 于 光速 时 ， 量 子 场 的 
方程 由 Klein-Gordon 方程 变 为 薛 定 廖 方程 ， 但 我 们 还 须 考 虑 标量 场 粒 子 之 间 的 引力 作用 ， 
因此 薛 定 廖 方程 带 有 一 项 引力 势 (8) 项 。 此 即 薛 定 谓 - 泊 松 方程 系统 (Schrödinger Poisson 
equations, SPE)， 有 具体 形式 为 : 


OW Rh? Los 

| | =0, 1 
ih EY jmd v + my =0 (10) 
V?$ó = —4rGp , p= m. (11) 


其 中 ， 式 (IU) whe 4 MARAE RY EYS 7 fe (NLSE). 

对 于 这 种 SPE 弱 引 力 系统 ，NLSE 方程 在 非 线 性 光学 中 已 经 存在 ， 因 此 不 用 担心 对 应 
关系 (参见 第 2 章 的 第 一 段 )， 重 要 的 是 寻找 泊 松 方程 的 光学 对 应 形式 。 这 里 我 们 介绍 介质 
的 热 光 学 效应 正好 具有 泊 松 方程 的 形式 ”。 

热 光 学 效应 指 光 场 中 的 介质 会 吸收 光子 能 量 ， 从 而 使 介质 温度 升 高 ， 继 而 导致 折射 率 改 
变 。 此 时 An = BAT, B 为 介质 的 热 光学 系数 ，AT — T — To, To 为 光 场 通过 介质 前 的 参 
考 温度 ， 为 固定 常数 。 实 验 上 达到 热学 稳 态 后 ， 介 质 中 的 热传导 满足 如 下 方程 : 


&V^T = k&V?AT = —ol , (12) 


其 中 ，A 为 介质 的 热传导 系数 ，a 为 介质 的 吸收 系数 。 需 要 注意 的 是 ， 式 (I2) 已 经 假设 
B2T1822 < OT, WAKE 2 轴 方 向 的 温度 改变 ， 宁 ? = a3 为 二 维 拉 普 拉 斯 算 子 。 为 
了 与 式 (E) 对 应 ， 将 式 (I2) 稳 态 的 热传导 方程 改写 为: 


V?An = = , I=. (13) 
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从 式 (3) 可 以 看 出 ， 稳 态 的 热传导 方程 即 泊 松 方程 ， 
) 与 (I)， 介 质 折 射 率 的 改变 An 类 比 于 引力 势 p» 


的 热 光 效 应 来 类 比 牛 顿 引 力 。 对 比 式 (ITI 

光 强 了 工 类 比 于 引力 系统 中 的 密度 po。 通过 合理 

的 SPE 天 文 对 象 。 

2.3 ”两 种 光学 仿真 原理 的 对 比 小 结 
EDINE 

行 比较 。 首 先 ， 仿 真 强 

力 SPE 系统 是 利 


E 


引力 场 利用 光子 流体 


ERR, $: 仿真 引力 系统 的 非 线 公 


了 强 引 力 场 和 弱 引 力 SPE 系统 的 光学 仿真 原理 
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因此 ， 我 们 可 以 很 方便 地 用 介质 


调节 这 些 光学 实验 参数 ， 我 们 就 可 以 仿真 具体 


后 ， 我 们 对 两 者 的 特点 进 
效 的 弯曲 时 空 度 规 ， 仿 真 弱 引 


IN. Al 
DSF 


构造 了 


j 热 光 效 应 来 类 比 牛 顿 引力 。 其 次 ， 这 两 种 实验 方案 所 需 的 光学 非 线性 
效应 并 不 相同 : 仿真 强 引 力 场 时 ， 折 射 率 变化 量 与 光 强 的 关系 为 An = 一 yw|? (b 为 激 》 


É 


脉冲 的 空间 分 布 轮廓 )， 对 于 y 而 言 这 是 一 种 局 部 的 非 线性 效应 ， 而 仿真 弱 引力 系统 需 
Qan = Pup, HI An = 六 RIwl?2dr， 这 是 一 种 非 局 部 的 非 线 怕 


x 


=) 
FH 


效应 ( 即 折射 率 改 变 


An Æ v 的 空间 积分 )， 此 处 细节 请 参照 42 小 节 。 并 且 ， 由 An 的 表达 式 可 以 看 出 ， 仿 真 
强 引力 场所 需 介质 为 自 散 焦 介质 ， 仿 真 SPE 系统 所 需 介质 为 自 聚 焦 介 质 。 

总 而 言 之 ， 根 据 待 仿真 的 天 文 对 象 的 不 同 特点 ， 研 究 者 采取 不 同 的 实验 方案 ， 保 证 天 文 
系统 和 光学 系统 满足 同样 的 数学 形式 ， 并 且 巧 妙 地 实现 这 两 种 系统 具体 参数 之 间 的 一 一 对 


3 仿真 强 引力 系统 的 光学 实验 


在 第 2 HH, RAINE 
动 类 比 测试 粒子 在 弯曲 时 空中 的 运动 。 在 本 
在 水 流 中 如 何 仿真 旋转 黑洞 ， 然 后 过 渡 到 光 
j 及 其 超 辐 射 的 研究 进展 。 在 水 流 的 流 场 中 ， 


y= 


了 光学 实验 仿真 强 引力 系统 的 原理 


， 即 用 声波 在 光子 流体 中 的 
章 中 ， 由 于 水 流 的 流 场 比较 直观 ， 我 们 首先 介 
学 系统 (光子 流体 )， 介 绍 光学 实验 仿真 旋转 黑 
不 仅 声波 而 且 表面 重力 波 (gravity wave) 的 传 


Sb AY wm 


£7; fitu 5 9| 73 P EE IR RC ER n Pe 8 


方程 相同 ， 即 可 以 用 重力 波 作为 流 场 


FP 的 扰动 ， 


来 类 比 测试 粒子 在 引力 场 中 的 运动 ， 从 而 描绘 


解 ， 下 面 我 们 就 用 重力 波 实验 来 介绍 仿真 原理 。 
作为 旋转 黑洞 的 Kerr 度 规 。 此 流 场 要 满足 2.1 小 节 所 提 


首先 要 构造 一 个 稳定 的 流 场 ， 
到 的 几 点 要 求 ， 即 无 条 的 势 流 ， 除 了 奇 点 之 
一 个 排水 口 ， 并 且 给 流体 一 定 角 动 量 ， 这 样 在 
被 称 为 “浴缸 漏 点 ”模型 吕 。 基 于 “浴缸 漏 点 
向 流速 w = 一 二 ， 角 向 流速 ve = —, 


BRR Dit ze Sh 
动 不 可 
红 移 面 与 视界 之 间 的 区 域 称 为 能 层 ， 声 波 或 
旋转 的 强烈 拖 点 效应 。 如 果 它 有 盏 没有 坠落 到 


F 视 界 )， 视 界 以 


A. BH 
根据 视界 和 能 层 的 定义 ， 流 场 的 径 向 流速 w 与 声速 c, 相同 之 处 是 模拟 黑洞 的 视界 ( 严 
内 流体 的 径 向 流速 均 超过 声速 ， 因 此 视界 内 流体 
能 通过 视界 向 外 传播 ， 流 体 的 速度 v SER c 相同 的 地 方 ， 


出 引力 场 度 规 。 水 流 中 的 重力 波 实 验 很 容易 理 


处 处 无 旋 ， 不 存在 湛 流 。 一 般 是 在 水 池 中 放置 
排水 口 附 近 就 形成 了 一 个 涡 旋 流 场 ， 这 个 模型 
”模型 ， 我 们 可 以 计算 流 场 的 流速 w， 其 中 径 
司 定 参数 ， 由 排水 速率 和 初始 角 动 量 决定 。 


的 声学 扰 
为 无 限 红 移 面 。 无 限 
E 力 波 等 扰动 在 能 层 区 域 能 感受 到 与 黑洞 一 起 
视界 内 ， 就 会 被 黑洞 散射 并 且 携 带 更 高 的 能 量 


46 天 文学 进展 42 卷 


逃逸 出 来 ， 这 就 是 超 辐 射 ” ”。 在 经 典 黑洞 研究 领域 ， 超 辐射 的 物理 本 质 是 提取 黑洞 旋转 
能 的 彭 罗斯 过 程 ” ， 因 此 超 辐射 有 时 又 被 称 作 彭 罗斯 超 辐 射 。 

Torres 等 人 ”利用 水 流 中 的 重力 波 来 测量 超 辐射 效应 ， 如 图 1 所 示 。 在 水 池 中 间 构 
造 出 稳 态 的 流 涡 流 场 ( 即 “ 浴 和 缸 漏 点 ”模型 所 类 比 的 旋转 黑洞 )， 然 后 从 左边 激发 出 重力 波 
扰动 ， 重 力 波 经 过 旋转 黑洞 时 产生 散射 波 ， 并 与 入 射 的 重力 波 发 生 干 涉 ， 生 成 明显 的 干涉 
图 样 。 若 经 过 黑洞 散射 后 ， 出 射 的 重力 波 强度 大 于 入 射流 ， 则 说 明 发 生 了 超 辐 射 效 应 。 实 
验 ” 表明: 当 入 射 重 力 波 的 角 动 量 与 中 心 流 涡 (类 比 旋转 的 黑洞 ) 同 向 时 ， 在 某 些 条 件 下 就 
可 以 产生 波幅 被 增强 的 出 射 波 。 


注 : 图 中 颜色 条 代表 高 度 ， 即 表面 重力 波 的 波幅 。a) 和 b) 分 别 是 不 同 频率 的 重力 波 (作用 等 同 于 声波 ) 干涉 图 


样 ， 由 此 干涉 图 样 可 以 计算 散射 后 出 射 波 的 幅度 。a) 中 频率 f — 3.45 Hz. b) 中 频率 f — 3.70 Hz. 


1 ”入 射 表面 重力 波 被 尝 涡 散射 后 ， 出 射 波 和 入 射 波形 成 的 干涉 图 样 ” 


3.1 ”光学 实验 仿真 旋转 黑洞 

与 上 述 水 流 相 比 ， 光 学 中 的 稳定 流 场 不 太 直 观 ， 但 是 依旧 可 以 参考 “浴缸 漏 点 模型 
测定 光子 流体 的 径 向 和 和 角 向 流速 。 参 照 2.1 节 ， 光 子 流体 的 速度 = Lo Vo. o 为 光 场 肪 
冲 包 络 的 相位 。 因 此 ， 在 光学 实验 中 为 了 构建 稳定 的 涡 旋 流 场 ( 即 光 涡 旋 )， 需 要 引入 相位 。 
该 实验 设想 最 初 由 Marino 四 提出 ， 并 由 Vocke 等 人 四 实现 ， 如 图 2 所 示 。 激 光束 在 经 过 相 
位 撼 模 后 ， 会 产生 多 级 衍射 光 ， 通 过 小 孔 将 一 级 衍射 光 (拓扑 荷 ! = 1 的 涡 旋 光束 ) 提取 出 
来 。 最 后 通过 分 束 镜 分 为 两 束 光 : 一 束 光 用 来 测量 光束 的 强度 分 布 ; 一 束 光 用 来 测量 光束 的 
空间 频谱 号 。 测 量 空间 频谱 时 需要 将 摄像 机 放置 在 透镜 的 焦 平面 处 ， 透 镜 的 作用 相当 于 对 
光 场 进行 二 维 傅 里 叶 变换 ， 在 焦 平 面 上 形成 的 像 就 是 透镜 前 光 场 的 空间 频谱 。 

光 涡 旋 的 波 前 呈现 螺旋 状 ， 携 带 相 位 因子 exp(i19)， 其 中 1 称 为 拓扑 荷 数 ，0 为 方位 角 ， 
即 光 涡 旋 绕 轴 旋 转 一 周 ， 相 位 改变 2x1。 因 为 其 特殊 的 拓扑 结构 ， 光 涡 旋 中 心 存在 相位 奇 
点 ， 这 就 造成 了 光 涡 旋 中 心 的 阴影 区 域 ( 光 场 强度 为 0)。 图 加 给 出 常见 的 三 种 产生 光 涡 旋 
的 方法 ， 第 一 种 方法 是 螺旋 相位 板 ， 它 是 一 个 沿 着 中 心 方位 角逐 渐变 厚 的 透明 板 ， 厚 度 与 方 


T 
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相机 5 
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相机 远 场 成 像 


4T 


EUN =i = 
Ar 


近 场 成 像 


i 


波长 532 nm 的 连续 激光 ， 通 过 快门 控制 激光 的 打开 时 间 ， 利 
一 2rVr/ro， 然 后 通过 小 孔 选 择 特定 级 次 的 涡 旋 光束 ， 输 入 到 
远 场 测量 光 强 分 布 及 空间 频谱 


仿真 旋转 黑洞 的 光 涡 旋 实验 装置 图 


A 小 孔 相位 掩 模 连续 激光 器 


注 : 实验 使 
TT p = 20 


分 别 在 近 场 和 


相位 掩 模 为 激光 束 产生 指定 的 相位 
I 介质 中 进行 演化 ， 在 介质 中 演化 后 


图 2 


位 角 成 正比 ， 类 似 一 个 旋转 台阶 ， 如 图 Ba) 所 示 。 第 二 种 方法 是 计算 得 到 相位 图 ( 见 图 
然后 加 载 到 空间 光 调 制 器 (spatial light modulator, SLM)。 第 三 种 方法 是 计算 全 息 法 ， 
是 Vocke 等 人 ”所 用 的 方法 。 通 过 计算 光学 算法 得 到 全 息 图 后 ， 可 以 将 全 息 
干 板 上 制 成 光栅 (如 图 区))， 也 可 称 为 相位 掩 模 。 全 息 图 也 可 以 加 载 到 空间 光 调 制 器。 


AO | 
相位 po 
b) 


强度 
输入 高 斯 光束 


输出 涡 旋 光束 


相位 


le} 


强度 


图 


È: a) 螺旋 相位 板 ; b) 相位 


加 载 到 SLM; c) 全 息 光 栅 。 


Bb)), 
也 正 


EST Ep 4 


图 3 ”制备 涡 旋光 束 常用 的 三 种 方法 
Vocke 等 人 ”通过 近 场 测量 光束 的 强度 分 布 ， 即 光子 流体 的 密度 p， 从 而 计算 出 流 场 每 
一 点 的 声速 d = 一 二 。 远 场 前 为 一 个 扫描 器 (图 中 的 小 孔 )， 它 沿 着 光束 的 直径 方向 扫 


可 以 测 


0 
量 每 一 点 处 的 空间 频谱 (Ks。,，K,)， 从 而 获得 流 场 每 一 点 的 速度 ， 计 算 公式 如 下 : 


(03 


通过 实验 测量 流 场 中 每 一 点 的 声速 及 流速 ， 可 以 确定 实验 仿真 的 弯 
黑洞 的 视界 及 能 层 ， 如 图 4 所 示 。 通 过 光子 流体 构造 出 等 效 的 时 空 结构 之 后 
的 问题 就 是 如 何 模拟 旋转 黑洞 的 超 辐射 过 程 ， 
KTT, 2022 年 Braidotti AM" 首次 报告 了 超 辐射 效应 的 测量 结果 ( 见 3.2 节 )。 


C 


c AK.) 


K,(r) 


noko noko 


, 


(14) 


和 时 空 结构 ， 识 别 出 旋 转 
， 一 个 很 自然 
该 研究 领域 在 Vocke 等 人 ”之 后 很 快 发 展 起 
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V/(105m-s?) 
V/(105m:s) 


r/mm 
a) 


YE: a) b) 分 别 为 黑洞 视界 及 能 层 的 计算 结果 与 实测 结果 ， 红 色 线 与 蓝 色 线 分 别 为 总 流速 v 和 径 向 流速 v.， 黑 
色 线 为 声速 cs， 前 后 两 个 交点 分 别 为 旋转 黑洞 的 视界 及 能 


4 ”光子 流体 中 的 旋转 黑洞 


p 


3.2 ”光学 实验 测量 旋转 黑洞 的 彭 罗斯 超 辐射 

旋转 黑洞 的 彭 罗斯 过 程 ， 已 在 前 文 作 了 简单 介绍 。Braidotti 等 人 ”将 一 束 强 度 较 弱 的 
SG (作为 光子 流体 中 的 声波 ) 与 泵 浦 光一 起 输入 介质 中 ， 如 图 5 fim. St SRE 
相互 作用 后 ， 产 生 了 一 束 新 的 闲 频 光 ， 并 且 被 困 在 能 层 内 ， 而 信号 光 成 功 地 逃 出 能 层 。 在 光 
学 系统 中 ， 信 号 光 汲取 泵 浦 光 能 量 的 过 程 也 被 称 为 光学 参量 放大 ， 但 这 一 过 程 对 相位 匹配 要 
求 非常 严格 。 具 体 而 言 ， 泵 浦 光 在 传播 过 程 中 不 断 激 发 每 一 处 信号 光 和 闲 频 光 ， 如 果 产 生 的 
着 号 光 与 前 面 已 经 存在 的 信号 光 相 位 不 一 致 ， 信 号 光 的 强度 就 可 能 会 减弱 。 只 有 新 产生 的 信 
号 光 与 前 面 的 信号 光 步 调 一 致 时 ， 才 能 最 终 增 强 信号 光 。 相 位 匹配 的 方法 有 很 多 种 ， 在 旋转 
黑洞 的 实验 中 ， 它 被 称 为 几何 诱导 相位 匹配 ， 即 闲 频 光 与 信 号 光 在 空间 上 分 离 ， 其 对 应 彭 
斯 过 程 ， 而 相位 匹配 条 件 也 与 产生 超 辐射 的 Zel'dovich-Misner 条 件 ” 一致 


o 


MARHA m 
答 入 信号 光 、 ‘int OIE 


性 介质 


图 5 ”仿真 彭 罗斯 超 辐 射 的 光学 实验 示意 图 
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最 后 ， 我 们 看 一 下 超 辐 射 的 测量 结果 ， 如 图 6a) Aras. KAP, hom q 分 别 为 泵 浦 光 、 
探测 光 和 闲 频 光 的 拓扑 荷 数 ( 详 见 3.1 78). Ry 为 反射 系数 ， 是 衡量 超 辐射 是 否 发 生 的 指标 ， 
Ry > 1 即 为 发 生 超 辐射 的 情况 。 图 6b) 显示 了 图 6a) 中 两 个 不 同 参 数 信号 光 与 闲 频 光 的 空 
间 分 布 ， 发 生 超 辐射 时 ， 信 号 光 与 闲 频 光明 显 分 离 ， 未 发 生 超 辐 射 时 ， 两 者 重 县 。 


P 


0 500 1000 1500 
r/um 
a) b) 


YE: a) 反射 系数 Ry 测量 结果 ， 泵 浦 光 的 拓扑 荷 数 1 = 2， 信 号 光 的 拓扑 荷 数 为 n， 闲 频 光 的 拓扑 荷 数 
q — 21 — n. b) 信号 光 与 采 频 光 的 空间 分 布 测 量 结果 ， 上 图 对 应 a) 中 q = 5 的 情况 ，Rvy < 1， 未 产生 超 辐 
射 ， 下 图 对 应 a) 中 q = 0 Wil, Ry > 1， 产 生 了 超 辐 射 。 


图 6 ” 超 辐射 测量 结果 
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4 仿真 弱 引 力 SPE 系统 的 光学 实验 


在 本 章 中 ， 我 们 介绍 两 个 仿真 弱 引 力 系 统 的 光学 实验 : Bekenstein 等 人 "利用 高 斯 光 
束 与 艾 里 光束 的 相互 作用 ， 模 拟 大 质量 天 体 产生 的 引力 透镜 等 效应 ，Roger A 通过 光 涡 
旋 在 介质 中 的 演化 ， 模 拟 旋转 “ 玻 色 星 ” ( 非 相 对 论 性 ) 的 演化 。 

4.1 ”光学 实验 仿真 引力 透镜 等 效应 
4.1.1 实验 方案 

为 了 模拟 大 质量 天 体 (比如 大 质量 恒星 、 星 系 和 星系 团 ) 的 引力 透镜 、 潮 汐 力 、 引 力 红 
移 和 蓝 移 等 效应 ，Bekenstein 等 人 ™ 设计 了 激光 实验 。 由 于 这 些 引 力 效应 属于 弱 引 力 场 现 
象 ， 无 需 模拟 完全 的 爱 因 斯 坦 场 方程 ， 因 此 他 们 采用 2.2 节 中 所 介绍 的 基于 热 光 效应 的 SPE 
系统 模拟 方式 。 

他 们 利用 波长 为 488 nm 的 连续 激光 ， 激 光 器 输出 的 高 斯 光束 经 过 分 束 器 后 分 为 两 束 


= 
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光 。 当 一 束 光 经 过 空间 光 调制 器 时 ， 产 生 一 种 特殊 的 光束 ， 叫 做 “加 速 光束 ”， 然 后 与 另 一 
束 高 强度 的 高 斯 光束 一 起 进入 介质 (光学 玻璃 SF11) 中 并 相互 作用 ， 以 此 来 模拟 各 种 引力 效 
应 ， 如 图 Ta) 所 示 。 


b) 


Hk: a) 实验 光路 图 ，b) 艾 里 光束 轨迹 示意 图 。 


图 7 ”仿真 引力 透镜 等 效应 的 光学 实验 


激光 器 输出 的 光束 是 在 谐振 腔 边界 条 件 下 的 特 解 ， 其 振幅 的 横向 分 布 为 高 斯 分 布 ， 因 此 
称 之 为 高 斯 光束 。 高 斯 光束 在 输入 介质 之 后 ， 由 于 热 光 效应 ( 详 见 2.2 节 )， 会 使 介质 中 心 折 
射 率 高 及 边缘 折射 率 低 ， 从 而 使 光 向 中 心 会 聚 ， 进 而 产生 一 种 自 聚 焦 效 应 。 另 一 方面 ， 高 斯 
光束 在 传播 过 程 中 会 有 衍射 效应 ， 宽 度 不 断 增 加 。 当 自 聚 焦 效 应 与 衍射 效应 完全 抵消 时 ， 高 
斯 光束 就 在 介质 中 稳定 地 传输 。 加 速 光 束 ， 又 称 艾 里 光束 ， 在 传播 时 具有 自 弯 曲 传输 ， 无 衍 
射 等 优点 。 由 于 艾 里 光束 这 些 优 点 ， 其 轨迹 及 结构 的 变化 可 以 很 好 地 反映 引力 作用 ， 不 会 受 
到 衍射 效应 的 影响 。 
4.1.2 仿真 实验 结果 

在 仿真 实验 中 ， 高 光 强 的 高 斯 光束 类 比 于 大 质量 天 体 ， 艾 里 光束 类 比 于 引力 场 中 的 测试 
光子 ; 前 者 所 产生 的 热效应 会 影响 艾 里 光束 的 轨迹 和 结构 ， 类 比 于 大 质量 天 体 所 产生 的 引力 
效应 。 艾 里 光束 在 介质 中 自 弯 曲 传输 ( 见 图 7b) 中 蓝 色 线 )， 由 于 高 斯 光束 ( 见 图 7b) 中 红色 
线 ) 的 影响 ， 艾 里 光束 的 弯曲 程度 减 小 ( 见 图 Tb) 中 绿色 线 )， 这 一 现象 就 类 比 于 引力 透镜 效 
应 。 随 着 激光 输出 能 量 的 提高 ， 我 们 可 以 观察 到 艾 里 光束 出 射 时 更 偏向 高 斯 光束 一 侧 ， 即 引 
力 透 镜 效应 的 强度 也 在 增加 ( 见 图 8a))。 

更 有 意思 的 是 ， 随 着 两 者 的 相互 作用 ， 艾 里 光束 的 结构 也 发 生变 化 ( 见 图 8b))。 并 且 随 
着 光 强 增 大 ， 结 构 的 改变 更 为 显著 ， 这 一 现象 类 比 于 大 质量 天 体 的 潮汐 作用 。 艾 里 光束 各 部 
分 与 高 斯 光束 的 距离 并 不 相同 ， 这 就 类 比 于 处 在 大 质量 天 体 周围 的 天 体 ， 其 各 部 分 所 受到 的 
大 质量 天 体 的 引力 场 并 不 相同 ， 从 而 受到 潮汐 作用 导致 变形 。 

除 此 之 外 ，Bekenstein 等 人 通过 控制 入 射 艾 里 光束 的 轨迹 ， 使 它 远离 或 者 靠近 高 斯 
光束 ， 并 且 测 量 出 射 时 艾 里 光束 波 矢 % 的 变化 。 他 们 发 现 ， 在 靠近 高 斯 光束 时 ，k; 会 增 
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HE: a) 引力 透镜 效应 ，b) 潮汐 效应 。 


图 8 ”引力 透镜 效应 与 潮汐 效应 加 


加 ， 反 之 则 减 小 。 在 引力 系统 中 ， 光 子 在 离开 大 质量 天 体 的 引力 势 阱 时 ， 会 发 4 


E 引力 红 


移 现 象 ， 因 此 光 的 频率 降低 ;反之 则 为 引力 蓝 移 ， 从 而 使 光 的 频率 升 高 。 光 学 实验 中 艾 


HGH ke 的 变化 与 光子 频率 的 变化 类 似 ， 折 射 率 An 的 变化 与 引力 势 $ 的 变化 类 似 ， 
2 = An(zs) — An(z1)， 因 此 k, 的 变化 就 是 仿真 引力 红 移 ( 蓝 移 ) 现象。 
4.2 ”光学 实验 仿真 旋转 “ 玻 色 星 ”演化 
在 天 文学 中 玻 色 星 全 部 由 玻 色 子 构成 。 恒 
力 收缩 ， 保 持 稳定 的 状态 。 与 此 不 同 ， 玻 色 星 是 一 种 玻 
观 量子 态 ， 由 于 测 不 准 原理 (或 者 说 波动 性 ) BEC 4 
章 )， 此 外 ， 玻 色 子 之 则 也 可 能 存在 排斥 性 自 相互 作 朋 
力 


Ed 
A 


Ed 
A 


色 凝聚 


Cy 


- 爱 因 斯 坦 


这 两 种 作用 可 以 平衡 玻 色 星 自 


H, 


即 


一 般 是 通过 核 聚 变 产生 的 辐射 压 来 阻挡 引 
体 (BEC), BFR 
HAPTER (或 称 作 衍射 效应 ， 见 第 2 


引 


H 


玻 色 星 可 以 旋转 ( 即 具 有 角 动 量 )，Roger 等 人 ”把 旋转 的 玻 色 星 进行 了 如 下 处 理 : 
星 整 体 为 一 个 涡 旋 ， 中 心 有 
论 极限 ， 即 SPE 系统 。 关 于 这 两 种 处 理 
会 提 到 。 
4.2.1 光学 实验 方案 

由 于 光学 实验 只 能 
为 时 间 轴 )， 因 此 NLSE 和 
其 温度 改变 量 (或 折射 率 改 变 


HAZ; FE Ne ARE ( 见 第 2 章 )。 对 于 二 维 泊 松 方程 
量 ， 对 应 于 引力 势 ) 的 形式 是 对 数 函数 ， 在 7 一 oo 时 不 趋 了 


里 ， 


玻 色 
个 相位 奇 点 。 另 外 ， 他 们 把 玻 色 星 当 作 是 EKG 系统 的 非 相 对 
(他 们 对 玻 色 星 的 “ 建 模 ”)， 我 们 在 4.2 节 结 尾 处 


研究 光 场 在 横 截 面 上 随时 间 的 变化 (2 轴 作 为 光束 传播 的 方向 ， 可 作 
系统 而 言 ， 
FO, 


因 


此 数值 解 泊 


松 方程 时 一 直 要 积分 到 边界 处 。 也 就 是 说 ， 在 实际 实验 中 ， 无 论 光 场 的 边界 有 


多 远 ， 边 界 条 件 对 于 引力 势 积分 的 影响 


都 不 可 忽略 。 因 此 与 3D 系统 相 比 ， 在 实际 工作 中 ， 


2D 系统 上 述 跟 边界 条 件 相 关 的 处 理 是 个 大 麻烦 (参考 文献 区 


的 2.4 节 )。 一 种 常规 的 处 理 方 


法 是 , HW 


白松 方程 中 添加 一 个 屏蔽 项 ， 


角 保 光束 光 强 分 布 范围 远 小 于 光 场 边界 ， 这 样 可 以 在 让 


(a1 - 5) wera) =-F 04) | (15) 
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从 而 使 积分 收敛 。 其 中 ，o 为 非 局 域 作 用 长 度 ， 一 般 可 以 取 o m D/2, D 为 介质 横 
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向 长 度 。 


Roger 等 人 号 在 研究 中 采取 了 上 述 边界 条 件 处 理 方法 。 他 们 提出 了 包含 边界 影响 的 “分 


布 式 热 损 失 模 型 ”(distributed loss model, DLM): 


OAT 


(poC) Ot 


K 
= KV? (AT) + al — AT , 


(16) 


其 中 ，o 的 含义 与 式 n) 中 的 相同 。 对 比 两 方程 ， 可 知 DLM 实际 为 带 屏蔽 项 泊 松 方程 的 
一 种 直观 理解 角度 ”。 在 实验 中 ， 只 需要 确保 输入 光束 的 空间 频谱 满足 条 件 : | 天 ,| > 1/0 


( 即 光 束 能 量 的 空间 分 布 不 能 过 于 广泛 )， 那 么 DLM Ær <o 时 就 是 二 维 泊 松 方程 的 解 。 
Roger 等 人 的 光学 实验 所 涉及 的 光学 原理 和 光学 技术 本 身 ， 与 3.1 节 中 所 介绍 的 光 涡 旋 
实验 ”相同 ， 区 别 在 于 仿真 的 对 象 以 及 仿真 对 应 方案 不 同 。 其 中 光 涡 旋 的 生成 技术 也 采用 


SIGH, FEARED AWW FEA FoR RES SFC. APO CHUA Pr. Ed 


文字 简要 介绍 了 光路 和 各 个 器 件 。 


中 的 说 明 


D. ye 
T b N 激光 器 
"m 2 JER a D 


PAS 相位 掩 模 


远 场 


铅 玻 璃 (SF6) 


ik: 实验 使 用 波长 532 nm 的 连续 激光 ， 利 用 衰减 器 控制 输入 光束 的 功率 ， 通 过 相位 掩 模 ( 详 见 图 回 ) 为 激光 束 


附加 相位 8 = e 产生 出 光 涡 旋 ， 输 入 到 玻璃 (SFO 是 一 种 热效应 比较 强 的 特种 玻璃 ) 中 进行 演化 ， 演 化 后 通过 分 


束 器 将 光 分 为 两 束 ， 使 用 近 场 测量 光 强 分 布 ， 远 场 测 量 空间 频谱 。 


9 ”仿真 “ 玻 色 星 ” 演 化 的 实验 装置 图 


4.2.2” 光 涡 旋 仿真 实验 结果 


介质 中 光 涡 旋 的 演化 如 图 中 所 示 ， 图 中 的 横 坐 标 可 作为 时 间 轴 。 图 10b) 为 空间 频谱 的 
演化 。 在 整个 实验 过 程 中 ， 光 束 始终 满足 4.2.1 节 所 述 的 条 件 (ok)? 六 1。 需 要 指出 的 是 ， 


在 实验 中 ， 光 涡 旋 在 介质 中 随时 间 的 演化 不 能 实时 监测 ， 只 能 监测 光束 从 介质 出 射 后 的 结 


(图 中 的 横 坐 标 ) 作为 衡量 传播 距离 z (即时 间 ) 的 指标 。 


这 个 仿真 实验 最 有 趣 的 结果 是 : 即使 在 光 强 很 大 时 (高 密度 时 )， 光 涡 旋 收缩 与 膨胀 交替 


果 ; 另外 ， 介 质 的 长 度 难 以 连续 变化 。Roger 等 人 通过 不 断 改 变 输入 光束 的 功率 ， 将 功率 p 


地 振 落 进行 ， 也 没有 由 于 自 肾 焦 效应 出 现 击 缩 ， 参 见 图 10a) 中 实 空间 光 强 分 布 T(r) 随时 间 
( 即 横 轴 ) 的 演化 。 这 里 的 热 光 自 聚 焦 效 应 ， 对 应 于 天 文 场景 中 在 自 引 力作 用 下 的 Jeans 914 
(参见 2.2 和 4.2.1 节 所 分 析 的 泊 松 方程 )。Roger FAUK HRA HANZ, EKA 


种 机 制 : 其 一 就 是 4.2 节 开 始 所 介绍 的 量子 压 ， 即 抗衡 自 引 力 、 保 持 玻 色 星 


稳定 的 


原理 ;这 
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归 一 化 强度 


输入 功率 /W 


ik: a) 光 涡 旋 强 度 分 布 的 演化 ，b) 光 涡 旋 空间 频谱 


b) 


10 JE (旋转 玻 色 星 ) 随时 间 的 演化 


个 机 制 与 角 动 量 无 关 。 其 二 是 光 涡 旋 本 身 的 拓扑 性 质 ， 他 们 认为 这 个 单 涡 旋 中 心 存 在 相位 奇 


mà. HES D BONG IN ACE. 

在 本 节 最 后 ， 我 们 对 Roger 等 人 的 仿 
题 。 该 实验 以 及 所 有 这 种 原理 的 光学 实验 ， 
于 在 常规 光学 实验 中 ， 光 子 在 介质 中 不 可 外 


真实 验 进行 如 下 评论 : (1) 他 们 的 仿真 对 象 有 问 
只 能 仿真 “2D 空间 +1D 时 间 ” 的 系统 。 这 是 由 
BAY “AE”, 而 光束 前 进 的 方向 (空间 轴 z) 对 


应 于 所 仿真 系统 的 时 间 轴 。Roger 等 人 的 光学 实验 声称 所 仿真 的 天 文 对 象 是 玻 色 星 ， 但 玻 色 


星 是 三 维 引 力 系统 ， 因 此 他 们 的 “仿真 对 应 关系 ”不 合适 。 相 应 的 二 维 引 力 的 天 文 对 象 可 以 


是 瞳 物质 纤维 状 结构 等 ， 我 们 将 在 5.2 节 简 单 介绍 。 其 实 ， 在 这 个 研究 团队 两 年 后 的 工作 中 
(BU 3.1 节 中 所 介绍 的 Vocke 等 人 的 仿真 实验 ”)， 他 们 明确 地 说 明了 只 能 用 光 涡 旋 实验 来 仿 


真 二 维 引力 对 象 (或 者 可 以 用 二 维 引力 来 


述 的 系统 )。(2) 具体 到 天 文 场景 中 的 玻 色 星 ， 是 


该 用 EKG 还 是 SPE 来 描述 ， 我 们 将 在 第 5 章 进行 讨论 。(3) RTI AAA 


个 解释 有 问题 ， 他 们 提出 的 第 一 个 机 制 量子 压 (或 称 衍 射 项 )， 并 不 能 始终 抗衡 自 引 力 ;， 关 于 


SPE 系统 的 这 一 点 ， 已 有 理论 分 析 ” ， 学 界 只 是 不 能 确定 这 种 夫 缩 ( 非 线性 演化 之 一 环 ) 的 


终点 产物 是 什么 。 他 们 提出 的 第 二 个 机 制 不 成 立 。 因 为 玻 色 星 作为 宏观 量子 态 ， 费 曼 ” 在 


1955 年 已 指出 来 ， 这 种 涡 旋 是 量子 化 的 ， 


不 会 是 单个 的 大 涡 旋 。 总 结 一 下 就 是 ，Roger 等 


54 天 文学 进展 42 卷 


人 提出 的 旋转 玻 色 星 模型 并 不 正确 ;特别 是 其 中 用 二 维 泊 松 方程 去 硬 凑 三 维 引力 系统 ， 导 致 
文章 中 有 些 地 方 行文 混乱 。 因 此 ， 这 个 实验 只 有 在 “单个 光 涡 旋 ” 这 个 层面 上 去 理解 ， 才 是 
正确 的 。 


5 MAMRE 


5.1 ” 强 引 力 场 的 仿真 实验 的 评述 

正如 我 们 在 第 1 章 所 介绍 ， 无 论 是 仿真 方案 的 理论 探讨 ， 还 是 实际 实验 的 开展 ， 最 先 
是 从 最 重要 的 问题 开始 ， 即 黑洞 的 霍金 辐射 ， 然 后 ， 扩 展 到 强 引 力 场 的 各 个 方面 。 由 于 强 引 
力 场 领 域 的 研究 涉及 黑洞 视界 附近 的 量子 效应 ， 以 及 量子 场 论 与 广义 相对 论 的 协调 ， 因 此 吸 
引 了 许多 物理 学 家 (包括 理论 物理 学 家 和 实验 物理 学 家 )。 仿 真实 验 手段 也 不 仅 局 限于 光学 ， 
如 第 1 章 所 提 及 的 量子 纠缠 实验 ， 以 及 张 保 成 等 人 所 介绍 的 冷 原 子 实验 ， 因 此 我 们 不 作 过 
多 总 结 和 展望 分 析 ， 请 读者 参见 文献 四 , 2 引 。 光 学 仿真 实验 参见 第 3 章 。 
5.2 SPE 光学 仿真 的 评述 和 一 个 展望 : 波动 暗物质 

关于 弱 引 力 场 SPE 系统 的 光学 仿真 实验 ， 我 们 在 第 4 章 做 了 介绍 。 而 且 ， 我 们 可 以 看 
到 这 个 领域 才刚 刚 兴 起 ( 近 10 年 内 )。 相 比 强 引力 场 的 效应 ，SPE 系统 与 传统 天 文 观测 有 更 
多 的 联系 (如 引力 透镜 、 引 力 红 移 所 关联 的 光学 图 像 和 光谱 观测 )， 因 此 我 们 在 此 做 进一步 
的 评论 。 特 别 是 ，4.2 节 所 介绍 的 实验 更 适用 于 当前 天 文 界 非常 感 兴趣 的 课题 : 波动 暗物质 
(waveDM)”。 下 面 我 们 专门 来 评述 和 展望 4.2 节 所 介绍 实验 的 天 文学 意义 。 

Roger 等 人 ™ 利用 NLSE 所 描述 的 光 涡 旋 ， 仿 真 旋 转 “ 玻 色 星 ”的 演化 。 对 于 真实 的 
天 文 场景 ， 玻 色 星 一 般 是 致密 的 、 被 自身 的 强 引 力 场所 束缚 的 天 体 ， 即 第 2 章 中 介绍 的 
Klein-Gordon Geon"", ， 因 此 应 该 用 EKG FEA 。 

对 于 SPE 版 的 Klein-Gordon Geon， 天 文学 领域 的 有 关 讨 论 非 常 活跃 ， 这 就 是 波动 暗 
物质 。 上 暗物质 是 当代 天 体 物理 与 粒子 物理 的 基础 性 问题 ， 目 前 的 标准 宇宙 学 模型 ， 即 “ 暗 能 
量 + 冷 暗物质 ”模型 (ACDM)， 可 以 解释 从 宇宙 微波 背景 到 星系 分 布 的 大 尺度 结构 ， 以 及 小 
到 1 Mpc 左右 空间 尺度 的 天 文 观 测 现象 。 然 而 在 星系 以 及 星系 周边 尺度 上 ， 还 一 直 存在 着 
若干 “小 尺度 挑战 ”“。 近 10 年 来 ， 波 动 暗物质 模型 引起 人 们 广泛 的 关注 ， 它 通常 指 质 量 
小 于 30 eV 的 暗物质 粒子 ， 取 星系 尺度 上 的 特征 速度 约 250 km/s， 则 其 德 布 罗 意 波长 可 以 
超过 粒子 的 平均 间距 ， 可 以 很 好 地 用 经 典 波 来 描述 ， 因 此 被 称 为 波动 暗物质 。WaveDM 与 
CDM 模型 相 比 ， 最 大 优势 在 于 有 可 能 解决 星系 尺度 上 “小 尺度 挑战 ” Schive 等 人 通过 
高 精度 的 数值 模拟 发 现 ， 在 相同 的 初始 条 件 下 ， 波 动 暗物质 框架 下 的 大 尺度 结构 与 CDM 的 
结果 相同 ， 但 在 星系 尺度 上 ， 几 乎 每 个 waveDM 学 的 中 心 都 存在 平头 结构 ， 中 心 的 密度 近 
似 恒定 ， 有 希望 解决 CDM 模型 中 的 “平头 - 尖 头 ”问题 ( 见 图 中 )。 

天 文中 的 暗物质 结构 ， 从 宇宙 大 尺度 结构 尺度 到 星系 尺度 ， 有 多 种 形态 。 在 10 Mpc 左 
右 尺度 的 宇宙 网 状 结构 里 : 有 了 暗物质 片 状 结构 ， 这 种 结构 对 近邻 其 他 物质 的 引力 作用 ， 可 
以 用 一 维 泊 松 方程 来 描述 ( 即 一 维 引力 系统 )”; 还 有 暗物质 纤维 状 结构 ”， 这 种 长 条 状 
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200 kpc 
Sr 


100 kpc 


a) 
YE: a) 星系 尺度 的 冷 暗物质 量 及 子 结构 可 ， b) 波动 暗物质 大 尺度 纤维 结构 、 晕 及 其 中 心 的 平头 结构 四 。 


11 ” 冷 瞳 物质 与 波动 暗物质 框架 下 的 暗物质 结构 (1 Mpc 及 以 下 尺度 ) 的 数值 模拟 


NS 


结构 对 近邻 物质 的 引力 作用 可 以 近似 为 无 限 长 的 质量 棒 ， 此 即 二 维 的 泊 松 方程 。 在 星系 团 
(1 Mpe 及 以 下 尺度 ) 以 及 星系 尺度 ， 则 分 布 有 大 大 小 小 的 暗物质 晕 和 子 巡 ， 以 及 各 种 子 结 
构 ”， 这 些 结构 则 一 般 需要 三 维 泊 松 方程 来 描述 。 

对 波动 暗物质 结构 的 相关 天 文 研究 ， 基 本 是 基于 NLSE 方程 的 数值 模拟 方法 ， 尚 未 开 
展 相 关 的 仿真 实验 。 但 数值 模拟 的 能 力 是 有 限 的。 首先 ， 对 于 一 些 复杂 的 动力 学 过 程 ， 如 波 
动 瞳 物质 的 波动 淇 流 ， 数 值 模拟 较 难 复 现 ， 再 例如 波 壕 缩 (wave collapse)， 己 超出 了 NLSE 
有 效 的 范围 。 其 次 ， 就 方程 本 身 而 言 ，SPE 的 求解 比 冷 暗物质 单纯 引力 泊 松 方程 求解 更 耗 
时 ， 因 此 导致 模拟 的 规模 (“粒子 数 ”) 和 精度 非常 有 限 ， 不 能 真实 描述 现实 场景 的 暗物质 
结构 。 
因此 ，4.2 节 的 光学 仿真 实验 为 我 们 提供 了 一 种 很 好 的 研究 思路 。 对 于 作为 三 维 引力 系 
统 的 波动 暗物质 (waveDM halo)， 如 我 们 在 4.2 节 最 后 提 到 的 ， 通 常 光学 实验 无 法 仿真 ， 
需要 另辟蹊径 的 实验 方案 。 对 于 二 维 引力 系统 的 波动 暗物质 纤维 (waveDM filament) 和 一 维 
引力 的 波动 暗物质 片 (waveDM sheet)， 则 通常 的 光学 实验 正 是 模拟 它们 的 动力 学 演化 的 好 
手段 。 特 别 是 ， 对 于 近年 来 天 文 观测 所 发 现 的 Mpc 量 级 的 旋转 暗物质 纤维 状 结构 ，4.2 
节 所 介绍 的 “单个 光 涡 旋 ” 实 验 正 是 它们 的 光学 对 应 体 。 


参考 文献 : 


] Barceló C, Liberati S, Visser M. Living Rev Relativ, 2011, 14: 1 

] 张 保 成 , RRE, 詹 明 生 . 中 国 科 学 : 物理 学 力学 天 文学 , 2014, 44: 879 

[3] Paredes A, Olivieri D N, Michinel H. Physica D: Nonlinear Phenomena, 2020, 403: 2301 
] 
] 
] 


[4] Hui L. ANNU REV ASTRON ASTR, 2021, 59: 247 
[5] Jafferis D, Zlokapa A, Lykken J D, et al. Nature, 2022, 612: 51 
[6] Gordon W. Ann Phys, 1923, 72: 421 


56 天 文学 进展 42 卷 


[7] Unruh W G. Phys Rev Lett, 1981, 46: 1351 

[8] Philbin T G, Kuklewicz C, Robertson S, et al. Science, 2008, 319: 1367 

[9] Rousseaux G, Mathis C, Massa P, et al. New J Phys, 2008, 10: 3015 

10] Lahav O, Itah A, Blumkin A, et al. Phys Rev Lett, 2010, 105: 2401 

11] Belgiorno F, Cacciatori S L, Clerici M, et al. Phys Rev Lett, 2010, 105: 0901 
12] Weinfurtner S, Tedford E W, Penrice M C J, et al. Phys Rev Lett, 2011, 106: 1302 
13] Faccio D. Contemporary Physics, 2012, 53: 97 

14] Steinhauer J. Nature Physics, 2015, 12: 929 

15] Müoz de Nova J R, Golubkov K, Kolobov V I, et al. Nature, 2019, 569: 688 
16] Torres T, Patrick S, Coutant A, et al. Nature Physics, 2017, 13: 833 

17] Vocke D, Maitland C, Prain A, et al. Optica, 2018, 5: 1099 

18] Prain A, Maitland C, Faccio D, et al. Phys Rev D, 2019, 100: 4037 

19] Solnyshkov D D, Leblanc C, Koniakhin S V, et al. Phy Rev B, 2019, 99: 4511 
20] Braidotti M C, Faccio D, Wright E M. Phys Rev Lett, 2020, 125: 3202 

21] Braidotti M C, Prizia R, Maitland C, et al. Phys Rev Lett, 2022, 128: 3901 
22] Barceló C. Nature Physics, 2019, 15: 210 

23] Almeida C R, Jacquet M J. European Physical Journal H, 2023, 48:15 

24] Marino F. Phys Rev A, 2008, 78: 3804 

25] FIG, 赵 娟 , 肖 振 中 , 等 . 激光 与 光电 子 学 进展 , 2022, 59: 1 

26] Kaup D J. Phys Rev, 1968, 172: 1331 

27] Penrose R. Riv Nuovo Cimento, 1969, 1: 252 

28] Zel'dovich Y B. Sov Phys JETP, 1972, 35: 1085 

29] Misner C. Phys Rev Lett, 1972, 28: 994 

30] Bekenstein R, Schley R, Mutzafi M, et al. Nature Physics, 2015, 11: 872 

31] Roger T, Maitland C, Wilson K, et al. Nature communication, 2016, 7: 1 

32] Liebling S L, Palenzuela C. Living Reviews in Relativity, 2023, 26: 1 

33] Feng S, Winful H G. Optics Letter, 2001, 26: 485 

34] Shen Y, Wang X, Xie Z, et al. Light: Science & Applications, 2019, 8: 90 

35] Goodman J. Introduction to Fourier Optics. Fourth ed. New York: Freeman Company Inc, 2017: 140 


36] Feynman R P. Application of quantum mechanics to liquid Helium. New York: Intersience Publishers Inc, 
1955 

37] Schützhold R, Unruh W G. Phys Rev D, 2002, 66: 4019 

38] Schunck F E, Mielke E W. Class Quantum Grav, 2003, 20: 301 

39] Schive H Y, Chiueh T, Broadhurst T. Nature Physics, 2014, 10: 496 

40] Bullock J S, Boylan-Kolchin M. ANNU REV ASTRON ASTR, 2017, 55: 343 

41] Chavanis P. Phys Rev D, 2016, 94: 3007 

42] Binney J, Tremaine S. Galactic Dynamics. Second ed. Princeton: Princeton University Press, 2008: 400 
43] Wang P, Libeskind N, Tempel E, et al. Nature Astronomy, 2021, 5: 839 


1 期 ERR, $: 仿真 引力 系统 的 非 线性 光学 实验 进展 57 


Developments of Nonlinear Optical Experiments Emulating 


Gravitational Systems in Astronomy 


CUI Jun-feng?, MENG Run-yu!'?, DONG Xiao-bo!?, ZHANG Xiao-shi?^^, 
LI Rui-xuan^?, BAI Jin-ming!? 


(1. Yunnan Observatories, Chinese Academy of Sciences, Kunming 650216, China; 2. University of 
Chinese Academy of Sciences, Beijing 100049, China; 3. School of Physics and Astronomy, Yunnan 
University, Kunming 650500, China; 4. Aerospace Information Research Institute, Chinese Academy of 
Sciences, Beijing 100094, China) 


Abstract: Nonlinear Schródinger equation (NLSE) is the common mathematical formula- 
tion of the laws governing many kinds of physical phenomena. In Optics, NLSE describe 
the spatio-temporal evolution of laser pulses in media. In Astronomy, it describe impor- 
tant astronomical objects and observed phenomena such as gravitational lensing, wave dark 
matter structures and substructures (e.g., cosmic filaments and haloes), etc. (the so-called 
Schródinger-Poisson equation systems, SPE), and even can be used to describe by analog 
the phenomena of strong gravitational fields (the Einstein Klein-Gordon systems, EKG). In 
recent years nonlinear optics has been developing rapidly, and optical experiments become a 
powerful method to solve or emulate NLSE. This article reviews the applications and devel- 
opments of the nonlinear optical experiments (the so-called “analogue gravity" experiments) 
for astronomical objects and phenomena: not only for several important observational SPE 
subjects such as gravitational redshift, tidal effects and (non-relativistic) rotating boson 
stars, but also for the hard-to-observe EKG phenomena such as black hole's Hawking ra- 
diation and Penrose superradiance. The analog/correspondence principle, optical design, 
experimental results and their astrophysical interpretations are presented. Finally we briefly 
outline the possibility of exploiting nonlinear optical experiments to emulate the dynamical 


evolution of various structures of wave dark matter. 


Key words: gravitational systems; wave dark matter; analogue gravity; experimental em- 


ulations; nonlinear optics; nonlinear Schródinger equation 


